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ABSTRACT 

We present results of Hubble Space Telescope (HST) observations of the 
radio galaxy 3C 244.1. The broadband F702W (R) and F555W (V) images 
(WFPC2/PC) show an elliptical galaxy and gaseous filaments and blobs sur- 
rounding it. In the narrow-band ramp filter, dominated by [O III]A5007, these 
filaments are bright and have the same morphology as the broad band images. To 
the south, the filaments have a cone-shaped structure and the radio jet is located 
at the center of this cone. To the north of the galaxy, the structure is found 
near the nucleus of the galaxy within its elliptical profile. From the photometry, 
the two brighter structures seem to be extended narrow line emission regions 
(ENLRs). The comparison with diagnostic line ratios shows that the observed 
emission is consistent with interactions between the expanding radio-jet and the 
local denser medium. 

Subject headings: AGN - Radio Galaxies - Jets 
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Association of Universities for Research in Astronomy, Inc, under NASA contract NAS 5-26555 
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1. Introduction 

The Hubble Space Telescope (HST) has been used to undertake a systematic survey of 
extragalactic radio sources. The SCR Snapshot Survey (de Koff et al. 1996; Martel et al. 
1998,1999 ; McCarthy et al. 1997, de Vries et al 1997) selected from the Revised Third 
Cambridge Catalogue of Radio Sources (SCR) (Bennett 1962a,b; Spinrad et al. 1985). These 
data, taken with the Wide Field Planetary Camera (WFPC2) using the F702W, F555W 
and ramp narrow-band filters, allow us to investigate the relationships between the radio 
and optical morphologies in a large sample of powerful radio galaxies over redhshifts of 
0<z< 1.5. 

Imaging of nearby Seyfert galaxies with HST, such as Mrk 3, Mrk 6, Mrk 573, NGC 
1068, NGC 2992 and NGC 4151, (Capetti et al. 1995a, 1995b, 1996, 1997; Winge et al. 1997; 
Allen et al. 1998 ,1999; Axon et al. 1998) has shown an intimate connection between the 
radio structure and the extended NLR (ENLR). These studies show that the interaction of 
the radio jet with the ISM is the main source of UV photons that ionize the ENLR in Seyfert 
galaxies. In previous work using HST data of SC 299, Feinstein et al. (1999) found that 
the the NE radio lobe lies within a shell-like structure of the ENLR, suggesting a physical 
connection between the jet and this ENLR. Evidence of this interaction comes from the 
values of the line ratio diagnostics and the [O II]A3727/[0 III]A5007 ratio, an estimator of 
the change of the ionization parameter f/, over the ENLR. We are extending this work by 
investigating whether this scenario is also applicable to other powerful radio sources. In this 
paper, we discuss the radio galaxy 3C 244.1 

3C 244.1 is a FRII radio source with two lobes with angular separation of 52" at po- 
sition angle (P.A.) of 168°. Both lobes have the same flux density (Fernini et al. (1997) 

log(L5GHz) ~ 22 W/Hz). Fernini et al. (1997) also detected the central component of the 
source that coincides with the optical identification of Spinrad et al. (1985), showing that 
the radio lobes are located asymmetrically with respect to this central compact source. 

At optical wavelengths, the host galaxy of 3C 244.1 was shown to have strong emission 
lines and a redshift oi z = 0.428 (Spinrad et al. 1985). From CCD images, McCarthy et al. 
(1995) note that the galaxy lies in a cluster, and their [O HI] image shows a high surface 
brightness fan of emission extending ~ 7" to the northwest, with a position angle similar to 
the radio structure. 

In this paper we discuss the morphology of SC 244.1 as shown by the high spatial 
resolution images of WFPC2, and we investigate the possible scenarios that give rise to the 
observed emission morphology. 
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2. Observations and data analysis 

The HST/WFPC2 observations of 3C 244.1 were taken as part of the 3CR imaging Snap- 
shot Survey (PI: Sparks) which was conducted in Cycles 4 to 8, in the F555W and F702W 
broad bands, and narrow emission-hne bands (de KofT et al. 1996; Martel et aL 1998,1999; 
McCarthy et al. 1997, etc). The 3C 244.1 images were obtained with the WFPC2/PC in 
April 1994 and May 1995 for the F702W filter and in July 1996 for the F555W filter. Narrow- 
band images with the WFPC2/WF2 and the FR680N ramp filter were taken in November 
1995. Table 1 shows the observation log. 

At a redshift oi z — 0.428, the galaxy is at a distance of 1.5 Gpc with a projected hnear 
scale of 5.1 kpc arcsec"^ (assuming Ho—65 km sec~^ Mpc~^ and Qo — 0.5). The WFPC2/PC 
scale is 0".0455 pixel"^ and for the WFPC2/WF2 the scale is 0".0996 pixel-^ and therefore, 
the physical scales for the images are 232 pc pixel"^ for the PC mode and 492 pc pixel"^ for 
the WF2 mode. 

The reduction procedure for the F702W filter data, which is close to Cousins R filter, is 
fully discussed in Martel et al. (1999). This reduction includes the standard WFPC2 pipehne 
processing followed by cosmic-ray removal. The data reduction was carried out using IRAF 
and the STSDAS package. At the redshift of 3C 244.1, the F702W filter includes the H7, 
[O III]A4363, He IIA4686, H/5, [O III]AA4959,5007 and [N I]A5199 emission hues. The F555W 
filter data, which is close to Johnson V, was reduced in the same manner and then registered 
and added to produce one final image with higher S/N. Due to the redshift, the F555W filter 
includes the flux from the emission hnes of [Ne V]A3426, [O II]A3727, [Ne III]A3869, RS and 
H7 (the latter considerably attenuated by the filter response) . In what follows we will define 
F702W as R, and F555W as V. The FR680N is a narrow-band ramp filter, which, given the 
redshift of the galaxy and its position on the CCD, produces an [O III]A5007 image. These 
images were reduced as described above and the final image was re-scaled to the resolution 
of the WFPC2/PC. 

The F702W and F555W filter data were fiux-cahbrated using values for the inverse sen- 
sitivities (PHOTFLAM) of 1.965 x 10"^^ ergs cm'^ A'^ DN-\ 3.49141 x 10"^^ cm'^ A'^ 
DN^^, and zero-points for the Vega system of Mphotzpt = —21.85 and Mphotzpt = —21.07, 
respectively. Non-rotated images were used for these flux calibrations to minimize interpo- 
lation errors. The FR680N filter was calibrated using the ramp filter calculator (Biretta 
et al, 1996). An approximate coordinate frame for the WFPC2 data is provided by the 
image header information, based on the HST guide stars. The accuracy for the coordinates 
is approximately 1" (Biretta et al. 1996). 

Images obtained with different filters were registered to a common reference frame using 
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cross-correlation of the brighter eUiptical structures. The sky background in each image was 
determined through the statistical analysis of a 3a clipping of the average background value 
in several regions in each of the images. 

3. Results 

The optical broadband images show an elliptical galaxy with blobs and filaments (Fig. 
1 and 2). Both the V and R images show these structures, and some are also present in 
the ramp image. We have named the most prominent of these structures as A,B,C,D and 
E (Fig. 3). To compare with radio data, we superimposed a 8.4 GHz radio map from the 
observations of Leahy (1999) in Fig. 1. The total extension of the radio structure is 51" 
(from lobe to lobe), larger that the size of these HST images. 

The object named A is bright and extended, and is located 2". 9 from the nucleus (P.A. 
125°), and has the same shape in all filters. The flux in the narrow-fllter implies a large 
emission of [O III]A5007 emission. The structure named B, is located 2". 6 away from the 
nucleus at P.A. 150°, and is similar in shape to structure A albeit fainter. Both of these 
structures (A and B) seem to be connected with the galaxy by very faint filaments. 

The object named C is located to the west at a separation of 1".4 from the nucleus at P.A. 
210°. A filament appears to connected it directly to the center of the galaxy. Superimposed 
over this line at 1" from the nucleus, there is another blob (C). The filament just described 
and those associated with blobs A and B, can be interpreted as the edges of a biconic structure 
with an opening angle of 85°. This kind of morphology has been found to be common in a 
variety of active galaxies ranging from nearby spirals (for example NGC 1068, Evans et al. 
1991; NGC 5728, Wilson et al. 1993) to powerful 3CR objects (e.g. 3G 299 Feinstein et al. 
1999). The southern radio-jet is located in the middle of this bi-conical structure. 

The structure D is located inside the elliptical profile of the galaxy, NE of the nucleus. 
After the subtraction of a r^/^ profile (using the IRAF/STSDAS task ellipse) from the image, 
a filamentary link emerges connecting D to the galaxy nucleus (see Fig. 1.). This structure 
seems not to be as large as the others and is very difficult to measure in the broad-band 
images because of the steep shape of the elliptical galaxy profile. This structure is also 
evident in the [O III]A5007 image, but as for the broad-band data, it is not possible to 
measure its fiux reliably. 

The structure named E, at 3". 3 from the nucleus at P.A. 110°, is not present in the 
narrow-band filter. Its profile is compact and steep, showing all the typical characteristics 
of a field galaxy with a redshift different from that of 3C 244.1 . 
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We have carried out photometry of the structures surrounding 3C 244.1 using the IRAF 
task polyphot (from the NOAO APPHOT package). This IRAF task computes the magnitude 
of an object in an image inside a polygonal aperture. The sky background is also computed 
and subtracted and the errors are calculated. To measure the bulk of the emission for each 
structure a polygon that closely matches the shape of the region was careful chosen. The 
same polygons were used on all images from the different filters to make the data comparable 
between these filters. However, for the V filter, where the signal is lower, some of these 
structures were difficult to measure above the background noise. 

Tabic 2 shows the result of the photometry using the polygonal measurement technique 
and the subsequent calibration. For the broad-bands filters, the data were caUbrated using 
the inverse sensitivity (PHOTFLAM) and the filter width (see Wide Field and Planetary 
Camera Instrument Handbook - Biretta, 1996). For the FR680N filter, the calibration 
was done using the IRAF/STSDAS synphot, and checked with the WFPC2 Exposure Time 
Calculator (ETC). The fiux observed at this filter, basically the flux of the [O III]A5007 line, 
was corrected by a small contamination (~ 12%) of the [O III]A4959 line (asumming that 
the flux of the [O III]A5007 hue is 2.88 times the [O III]A4959 flux). This correction was 
calculated by computing the central wavelength at the position on the CCD, and assuming 
that the filter is Gaussian, centered at that wavelength with a FWHM of 1.3% (value obtained 
from the ETC) of this central wavelength. Columns 2, 3 and 4 of Table 2 are the results of 
these calibrations. 

From Fig. 1, structures A and B are very close the jet flow, so the flrst attempt was 
to test if these objects were optical synchrotron knots on the radio jet. So we flt a power 
law to fiux from the broad band filters F555W and F702W. If the main source of fiux is 
synchrotron radiation, the fiux must follow a power law with an exponent in the range -0.5 
to -1.5 . To make the test as more accurately as we can, we used the task calcphot- synphot 
(from the SYNPHOT/STSDAS IRAF package) to integrate the power law under a better 
simulation of each filter. We found that the observed flux is not a power law and the best fit 
is a flat distribution with a null slope. So, we conclude that the observed flux is inconsistent 
with power-law distributions of optical synchrotron radiation. 

From the the broad band fllters, we can estimated the (V-I) colors of the elliptical host 
galaxy, A and B regions. Colors of the other structures can not be measured unless we have 
a better modeling of the underline galaxy flux, for which we need deeper observations and 
better spatial resolution. Objects A and B are far away from the elliptical proflle and are 
clearly not contaminated. For the elliptical galaxy we obtain (V-I) = 1.2 ±0.3 which is close 
to the value predicted by Ronnback et al. (1996) for an elliptical galaxy at the redshift of 
3C 244.1 (z=0.428) meanwhile the objects A and B, seems to be bluer with values -0.2 ± 
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0.2 and -0.6 ± 0.5 respectively. 

Both structures A and B, given the large flux in the narrow ramp fllter, are obviously 
at the 3C 244.1 redshift. Also both show the same morphology shape as the broad-band 
images. Since the R filter data also includes the fiux of H7, [O III]A4363, He IIA4686, H/3, 
[O III] AA4959, 5007 and [N I]A5199 emission lines, the total flux measured in this broad band 
must be larger than the ramp fllter which only includes [O III]A5007. The fact that they 
are so similar can be explained if we assume than the flux originates in a pure emission line 
region, without any signiflcant continuum contribution. 

As the narrow-band fllter data shows, there is a large amount of flux in the emission 
line [O III]A5007. Thus it is very reasonable to suggest that the structures A and B are 
Extended Narrow Line Regions (ENLR). The ENLR have been found to be associated with 
a wide variety of AGN, from Seyfert 2 to radio galaxies. 

Column 5 (Table 2) shows the flux calibration of the broad-band R fllter assuming 
that all the flux is from a source that has a spectrum dominated by bright emission lines. 
The shape of the F702W filter is practically fiat (See Wide Field and Planetary Camera 2 
Instrument Handbook, page 217), and so it is easy to compute the calibration of the fiux 
of all the lines by calculating the total fiux for only one emission line (that accounts for 
all the fine emission flux). To make this calculation, the ETC was used and the quantum 
efficiencies (system -|- fllter -|- CCD, hereafter QT) for each wavelength of the conspicuous 
emission lines included in the R fllter band were computed. Table 3 shows the results of this 
calculation and conflrm that all these emission lines are practically attenuated by similar 
amounts. Therefore we used H/5 as an equivalent line (i.e. a line that has all the flux of all 
the lines included in the fllter) for the calibration. The error of this procedure is extremely 
small. 



4. Discussion 

A key question is to identify the mechanism responsible for the ionization of these 
structures and fllaments. The most likely mechanisms are : photoionization by the AGN, 
photoionization by the AGN with matter-bounded clouds and shock-ionization by the radio 
jet. We discuss each of these mechanism in detail. 
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4.1. Photoionization by the AGN 

This mechanism assumes that the UV flux produced by the AGN's nuclear source pho- 
toionizes the line emission region. The physical status of the gas at any place can be described 
with the ionizing parameter (U), which considers the dilution of the UV photons as the dis- 
tance from the nucleus increases. For regions A and B, it is very unhkely that the AGN is 
the only source of photons because the large distances. To account for this scenario, some 
reference numbers of the total photon budget can be calculated. 

Rawlings et al. (1989) performed spectroscopy of the nuclear zone of 3G 244.1 and show 
a relation of [O 111]A5007 /H/? = 4.5. Using the results of Ferland and Netzer (1983), which 
assume that the main source of photoionization continuum is a power law {fi, oc and 
that the gas has a density of Nh — 10^, they calculate several line emission ratios. Using 
these results, it is straightforward to compute the ionization parameter ([/). This leads 
U = 3.2 X 10^"' for solar abundances and U = 6.3 x 10~^ for subsolar abundances. This 
spectrum covers only light coming from the nuclear zone that is observed in the narrow band 
filter (regions A and B are not included in this spectrum), so it is possible to calculate the 
total UV photon emission of the nucleus. Because Q = UAnr'^NeC and r can be estimated 
from the WFPG2 emission line image, this gives Q — 4.5 x 10^^ photons sec~^ for solar 
abundances and Q = 9 x 10^^ photons sec~^ for sub-solar abundances. Due to the large 
distances, these two values (for each this different abundances) are not enough to ionize the 
blobs A and B. Therefore, we conclude that direct photoionization by the nucleus is unlikely 
to be the dominant mechanism responsible for the observed ENLR emission. 



4.2. Photoionization by the AGN with matter-bounded clouds 

This idea was developed by Binette et al. (1996) in order to explain the discrepant 
high ionization line ratios and high electron temperatures observed in many active galaxies. 
These models are based on a parameter Am/i which is the solid angle covered by the "matter- 
bounded" component relative to that covered by the "ionization-bounded" component. In 
another paper, Binette et al. (1997) refined their models. Taking into account the possibility 
of having a larger U parameter, they derived three models with values oiU = 0.5, 0.05, 0.02 
and a density He = 1000 cm~'^ for the "matter bounded" clouds, which are exposed to 
ionization radiation from the AGN. 

As we only have one line and the integrated total fiux of several lines from the R filter, 
we computed the fiux ratio Co = R/[0 111]A5007 noting that the bulk of the R filter fiux 
is the sum of H7, [O 1I1]A4363, He 11A4686, H/?, [O I1I]AA4959,5007 and [N I]A5199 fiuxes. 
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Note that the contribution from the stellar host galaxy is not significant at the location of 
the blobs. The line ratio can be derived from Table 2 and are Cij = 2.6 for structure A and 
Cb = 2.0 for structure B. 

We computed Co (the ratio of hues include in the R filter to [O III]A5007) for the 
results of mixed models of the papers mentioned above. Table 4 shows these results of Co as 

a function of the Am/i parameter for the model of Bincttc ct al. (1996) and models H,M,L 
of (Binette et al. 1997). Models M and L seem to fit the observed Co for structures A and 
B, for Am/i ~ 0.2 and 0.34 for each object. But both of them have the same problem that 
we have discuss for the first case of photoionization by the AGN, U = 0.05 for Model M and 
U — 0.02 for model L, values that imply a very high Q {Q = cnJJdP) due to the distances 
to the galaxy nucleus. This calculation makes the mixed models very unlikely. 

4.3. Shocks 

Fig 1 shows that the ENLR structures A and B are located close to the edges of the 
radio-jet and it is possible that there is an interaction between the jet and these regions. It 
is important to study whether there is a physical relationship between the radio-jet and the 
optical hue emission. Taylor et al. (1992) proposed that fast bowshocks resulting from the 
interaction of the radio jet and the ISM were the source of ionizing photons of the emission- 
line gas in a number of sources. Capetti et al. (1995a, 1995b, 1996, 1997) and Winge et 
al. (1997) were the first to show that this mechanism best explains the optical emission in 
the NLR in nearby Seyfert galaxies (Mrk 3, Mrk 6, Mrk 573, NGC 1068, NGC 4151 and 
NGC 7319). Recently, this work has been confirmed by Aoki et al. (1999) and Kukula et 
al. (1999). In the case of the powerful radio galaxies of the 3CR catalogue, Feinstein et 
al. (1999) showed that this interaction also occurs in 3C 299, where clearly the NE radio- 
jet and the ENLR have similar morphologies, and where there is further evidence of this 
interaction from the values of the different emission-line ratios and the evolution of the line 
ratio [O II]A3727/[0 III]A5007, as an estimator of the changes of U, over the region. 

Dopita & Sutherland (1995a,b) have modeled in detail the ionization of the ENLR due 
to shocks. In one scenario, which has been shown to work for Seyfert galaxies, the radio jet 
interacts with the local interstellar medium and shocks the gas. In this scenario, the hot 
post-shock plasma gas produces photons that can diffuse upstream and downstream of the 
jet. Photons diffusing upstream can encounter the preshocked gas and produce an extensive 
precursor HII region, while those traveling downstream will influence the ionization and 
temperature structure of the recombination of the shock. 
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To check the vahdity of this interpretation, we use the hne diagnostic test of Dopita 
& Sutherland (1995b). We apphed the same procedure as the case above and use the flux 
ratio Cu = R/{0 III]A5007. On the other hand, we can compute the value of this ratio for 
emission lines arising in a shocked ENLR. Using Table 1 of Dopita & Sutherland (1995b) we 
calculated the total integrated flux for the lines involved and the ratio to [O III]A5007 flux 
as : 

[O7//]A5007 QT[oiii]xmm 

where i represents each of the lines mentioned above and their respective flux. 

Table 5 shows the results of this computation, for shock speeds of 100, 150, 200 and 500 
km sec~^, with and without a precursor H II region. Prom this table, the values observed 
for structures A [Co — 2.6) and B [Co — 2.0) are easily explained with shock models of 
~ 200 km sec~^. We conclude that shock interactions between the expanding jet and the 
local denser medium are responsible for the observed ENLR emission in 3C 244.1 



5. Conclusions 

We have shown that radio galaxy 3C 244.1 has a filamentary structure and some blobs 
bright in emission line ([O III]A5007). To the south, this structure seems to be an emission- 
line cone with an opening angle of 85° (from RA. 125° to RA 210°). The radio-jet is located 

at the center of this cone (R.A. 168°). The two brighter blobs (named A and B) arc likely 
part of the ENLR (by comparing the flux between broad-band and narrow-band), similar 
to those associated with AGNs. These structures arc larger and located far away from the 
nucleus (2". 9). To the north more of these structures were found, but near and associated 
with the nucleus. 

We tested several scenarios to explain the source of energy for the emission lines observed 
and we found that the direct photoionization by the AGN can only be possible if an unlikely 
large amount of UV photons is provided by AGN. The same is true in more sophisticated 
mixed medium models (Binettc ct al. 1996,1997). The observations seem to fit well the 
behavior of a model material shocked by the radio jet (Dopita et al. 1996b). A shock 
velocity of 200 km sec~^ can explain the ratio of the fines included in the i2-band filter to 
the [O III]A5007 fine-emission. 

Therefore, from the morphology (location of the radio-jet) and these physical arguments 
(the flux measured is consistent with emission from shocked gas) , we conclude that the ENLR 
structures of 3C 244.1 are the result of the interaction of the radio jet with ISM gas. 
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Table 1. Log of Observations 



Filter Name 


Emission Line 


Exp. Time 


Date 






sees 




F702W 


Broad band^ 


300 


1994 Apr 17 


F702W 


Broad band^ 


300 


1995 May 28 


FR680N 


[0 I1I]A5007 


300 


1995 Nov 13 


FR680N 


[0 111]A5007 


300 


1995 Nov 13 


F555W 


Broad band^ 


300 


1996 Sep 7 


F555W 


Broad band^ 


300 


1996 Sep 7 



^Duc to the filter width, the measured flux includes emi- 
sion from from several lines, see text 



Table 2. Total Flux^ from the Structures 



Region 


F702W 


F555W 


[0 III]A5007 


F702W^ 


A 


17.26 ± 1.37 


20.57 ± 2.69 


9.36 ± 1.24 


24.12 ± 1.91 


B 


5.25 ± 0.82 


5.04 ± 1.62 


3.58 ± 0.77 


7.34 ± 1.14 


C 


2.6 ± 0.59 






3.62 ± 0.83 


C 


2.3 ± 0.52 






3.24 ± 0.74 



^Flux units are 10 ergs cm ^ sec ^. 

^The F702W filter flux is calibrated assuming that the flux distri- 
bution is completely dominated by emission lines. See text. 



Note. — Measurements with values less than three times the error 
are not reported. 
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Table 3. Quantum efficiency (system + filter + CCD) for F702W at tfie line wavelength 



Line 


QT 


H7 


0.133 


[0 III]A4363 


0.139 


[He II]A4686 


0.132 


H/3 


0.124 


[0 III]A4959 


0.113 


[0 III]A5007 


0.112 


[N I]A5199 


0.092 



Table 4. Photoinization by an AGN with matter-bounded clouds 



Am/i A*^ H M L 

U = 0.04 U = 0.5 U = 0.05 U = 0.02 



,04 


1.79 


1.83 


4.95 


6.77 


14 


1.74 


1.76 


2.99 


3.19 


,24 


1.71 


1.71 


2.46 


2.50 


34 


1.68 


1.67 


2.21 


2.21 


,44 


1.66 


1.65 


2.06 


2.05 


54 


1.65 


1.63 


1.97 


1.95 


64 


1.64 


1.61 


1.90 


1.88 


,74 


1.63 


1.60 


1.85 


1.83 


,84 


1.62 


1.59 


1.81 


1.79 


,94 


1.61 


1.58 


1.78 


1.76 


,04 


1.60 


1.57 


1.76 


1.73 



^Model from Binette et al. (1996), no special 
name assigned in the original paper. For other mod- 
els of Binette et al. (1997), names are as in the 
original paper. 
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Fig. 1.— 3C 244.1, i^-band image taken with WFPC2/F702W filter (grey scale), the contours 
are from the radio map (Leahy, 1999) 



Table 5. Shock Models 



Shock Velocity 


Cd 


km sec~^ 




150 


2.29 


200 


2.67 


200 + precursor 


2.34 


300 


3.34 


300 + precursor 


1.67 


500 


2.71 


500 + precursor 


1.47 



Fig. 2.— 3C 244.1, top-left : i?-band image (F702W), top-right : F-band image (F555W), 
bottom-left : i?-band image with the elliptical profile substracted (note the plume of emission 
to the northwest), bottom-right : [O III]A5007 (re-scaled from the WFPC2/WF to the 
resolution of the WFPC2/PC) 



Fig. 3. — 3C 244.1, i?-band image with the identifications of the structures. The straight 
line indicates the direction of the radio jet. 
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